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Disclaimer 



Although this program was tested by its author prior to submission, 
no warranty, express or implied, is made by the author, 1620 USERS 
Group, or IBM as to the accuracy and functioning of the program and 
related program material and no responsibility is assumed by the 
author, 1620 USERS Group, or IBM in connection therewith. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter or pencil, do not use ink) 

Program No. Date 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No. 

it does? 

Comment _________ 

2. Does the program do what the abstract says ? Yes No 

Commen t 

3. Is the Description clear, understandable, and adequate? Yes No 

Commen t 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No 

Comment 

Are the Sense Switch options adequately described (if applicable)? Yes No 

Are the mnemonic labels identified or sufficiently understandable? Yes No 

Comment 

5. Does the source program compile satisfactorily (if applicable)? Yes N*o 

Comment 

^ 6. Does the object program run satisfactorily? Yes No 

C omment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No 

Comment 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No 

Comment 

9. Were all necessary parts of the program received? Yes No 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt Company 
Data Corporation 

7500 Old Xenia Pike Address 
Dayton, Ohio 45432 



User Group Code_ 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED 
TO PARTICIPATE IN THIS EVALUATION. 
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40-40 Correlation 



by 



Lew D. Harkins 

Computer Research Center 
Louisiana State University 
Baton Rouge, Louisiana 

November 10, 1961 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 







LOUISIANA STATE UNIVERSITY COMPUTER RESEARCH CENTER 

40 - 40 
Correlation Program 



Lew D. Harkins 



The 40 - 40 Correlation Program gives means, standard deviations, 
and correlation coefficients for all possible combinations of variables 
up to 40. Variables are entered in fixed point form and may be between 
2_and_7 digits in size. Sums of squares and crossproducts are accumulated 
and corrected (see Remarks) in fixed point arithmetic before being floated 
for the final calculations, thus reducing truncation errors inherent in 
floating point summations. The fixed point sums of squares and cross- 
products are limited to 16 digits; sums are limited to 10 digits, and n 
is limited to 5 digits. 



Input Format 

One observation on all variables to be correlated is entered on 
from one to four data cards, the size of the variables determining how 
many data cards will be needed to constitute one observation. Variables 
may be placed in any desired order on the cards, as their location is 
provided for the program's use by means of the parameter cards. A 
variable must be at least 2 .digits in size and not larger than 7 digits. 
It is not necessary to flag (x punch) the high order positions of the 
variables on the data cards. Negative numbers are denoted by a flag 
(x punch) over the low order polTition; "positive numbers, by the absence 
of a flag. 



Parameter Cards 

Three parameter cards are needed for each problem to be processed. 
All fields on these three cards are t wo digit fiel ds, and must_be flagged 
in the high order if used. All fields ~snus-t be used on parameter card 1, 
and only as many fields as there are variables must be used on each of 
cards 2 and 3. Unused columns on these three cards may be left blank. 



-1- 



Col. 


1 


- 2 


Col. 


3 


- 4 


Col. 


5 


- 6 


Col. 


7 


— 8 


Col. 


9 


- 10 


Col. 


11 


- 80 



Parameter Card 1 
Total number of variables 

Number of variables In first card of an observation 
Number of variables in second card of an observation • Zero if 
Number of variables in third card of an observation * fewer c 
Number of variables in fourth card of an observation / needed 
Blank observe 



Parameter Cards 2 & 3 

These two parameter cards specify the location of the variables on 
the input data cards. The terms "high order position" and "low order 
position" refer to column number on the data cards, without regard to 
card number, as the number of variables to be taken from each card is 
determined by parameter card 1. 









Parameter 


Card 2 








Parameter Card 3 






Col. 


1- 


2 


High order 


position 


of 


var. 


1 


Low order position of 


var. 


1 


Col. 


3- 


4 


High order 


position 


of 


var. 


2 


Low order position of 


var. 


2 


Col. 


5- 


6 


High order 


position 


of 


var. 


3 


Low order position of 


var. 


3 


Col. 


75- 


76 


High order 


position 


of 


var. 


38 


Low order position of 


var. 


38 


Col. 


77- 


78 


High order 


position 


of 


var. 


39 


Low order position of 


var. 


39 


Col. 


79- 


80 


High order 


position 


of 


var. 


40 


Low order position of 


var. 


40 



Input Procedure 

The 40 - 40 Correlation Program tests the last card indicator on the 
1620 to determine if all data for a problem has been read, therefore, 
separate problems cannot be stacked in the read hopper of the 1622 all at 
one time. If parameter cards and data for the first problems are stacked 
on top of the program deck it is only necessary to depress the start key 
after the program has been loaded, and the halt code (48) appears in the 
op register, to begin processing the first problem. For additional 
problems, it is only necessary to place parameter cards and data cards in 
the read hopper and depress the reader start key. 

If data cards are not stacked on top of the program deck, then the 
reset key must be depressed before depressing the start key after the 
program has loaded and the halt code appears in the op register. The 
program is self-initializing after each problem and need not be reloaded 
to process additional problems. The order of input for each problem is: 
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1) Parameter Card 1 

2) Parameter Card 2 

3) Parameter Card 3 

4) Data Cards 



Console Checks 



All tab stops on the typewriter should be cleared and only 
margins set. 

Program Switches 2, 3, and 4 are not used; program switch 1 

controls output (see "output") . 

Clear memory before loading the program. 

Input/output check switch should be set to stop. 

Parity check switch should be set to stop. 

Overflow check switch should, be set to stop. 

Digit located at core storage position 401 determines 

course of action in case of characteristic overflow or 

underflow, (see 1620 SPS Subroutines) , and may be changed 

after the program is loaded. Error codes are those of the 

SPS subroutines . 



Output 

Output is under control of program switch #1; with #1 on output 
is on cards; with #1 off output is" by typewriter . — — ~ 

Since the program treats all data as whole numbers, the location 
of the decimal in the output must be adjusted whenever the input data 
contain decimals. This is done as follows: 

means - output has one more decimal place than input data 
standard deviations"-" convert from floating point output to 
fixed point} then move the decimal point to the left the 
number of places in the input data. For example, if a 
standard deviation output is 5238456219 and the input for 
this variable contained 3 decimal places, then the standard 
deviation is .038456219^ correlati on coefficien t - no adjust- 
ment of decimal required. -™— — 
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REMARKS: 



Corrected as used here means corrected for the mean 

_ .2 

Ex -& ; Lxy - 

n n 

The following calculations are used to arrive at the means, 
standard deviations, and correlation coefficients: 




Floating point notation used is as follows: 



5018640000 = .1864, 5113700000 = 1.37, 4910000000 = .01, etc. 
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EXAMPLE PROBLEM 



5 Variables 
13 Observations 
(Raw Data) 
Observation X, X X 



1 


25.6 


2.9 


3. 


254. 


26.41 


2 


78.1 


6.2 


7. 


780. 


27.81 


3 


43.5 


5.0 




436. 


32.90 



Input Format 



. . 1st Data Card Last Data Card 

Col. 1-7 Blank Blank 

8-10 25,6 035 

11-12 29 Q8 

13-15 003, 008, 

16-18 254. 036. 

19-22 2641 4981 

23-80 B'lank Blank 



Parameter Cards 







1st 


2nd 


3rd 


Col. 


1- 2 


05 


08 


To 


Col. 


3- 4 


05 


11 


12 


Col. 


5- 6 


00 


13 


15 


Col. 


7- 8 


00 


16 


Is 


Col. 


9-10 


00 


19 


22 


Col. 


11-80 


Blank 


Blank 


Blank 



Output from Typewriter 



VR MEAN STD DEV 

01 41.48 5331493954 

02 437 5227424137 

03 40 5132914028 

04 4148 5331467761 

05 41042 5424G64444 

VR VR COR COEFP 

01 02 5096544667 

01 03 5060436826 

01 04 5099999007 

01 05 5011697431 

02 03 5050223114 
02 04 5096516763 

02 05 5016535123 

03 04 5060545014 

03 05 5013230260- 

04 05 5011827088 

Output Adjusted for Decimals 

VR MEAN STD DEV 

01 41.48 31.498954 

02 4.37 2.7424137 



05 41.042 24.864444 

VR VR COR COEFF 

01 02 .96544667 

01 03 .60436826 



04 05 .11327088 



Operating Instructions 



Clear memory. 



a. Set all check switches to program. 

b. Depress Instant Stop (SCE) , Reset, Insert 

c. Type: 26 00008 00009 

d. Depress Release, Start 

e. Allow 2 seconds, then depress Instant Stop (SCE) 

f. Set all check switches to stop. 



Load Program. 



a. Depress Reset (on console) 

b. Place program deck in read hopper 

c. Depress Load Key on 1622 unit 

d. To process last card, depress Reader Start Key. 

After program is loaded, the halt code (48) will appear in the operation 
register. Depress Reset, Start. (Reset Key must be depressed before the 
Start Key as the program interrogates the last card indicator.) 

If typewriter output is to be given, all tab stops should be cleared 
from the typewriter, and only margins set. 

Program switches should be set according to the program specifications. 

The program is now ready to accept parameter cards and data cards . 

The program is self- initializing between problems and need not be 
reloaded or restarted to process subsequent problems; however, the last 
data card in each set must be processed before the next set is placed in 
the read hopper (due to use of last card indicator). 
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1 INIT TFM *+30.6879Z 02.178 16 02208 06879 FLAGS 

2 AM *+18.16Z 02190 11 02208 00016 FLAGS 

3 TF 0.B1GB+14Z 02.202 26 00000 04741 

4 CM *-6.19999Z 02214 14 02208 19999 FLAGS 

_„? . BN1 * -36. 1200 Z 02226 47 02190 oi?nn _ 

6 TFM *+30.6479.9Z 02238 16 02268 06479 FLAGS 

7 AM *+18.10.9Z 02250 11 02268 00010 FLAGS 

8 TF SUMS.BIGB+BZ 02262 26 06489 04735 

9 CM *-6.6879,9Z 02274 14 02268 06879 FLAGS 
10 BNI *-36.1200Z 022B6 47 02250 01200 

J_L II* !±V. .?_ 0229S 16 04 976 0000 ..FLAGS 

12 RN R1-1.500Z 02310 36 04987 00500 

13 RN R1+9«500Z 02322 36 04997 00500 

14 RN R1+89.500Z 02334 36 05077 00500 

15 READ TFM RX1+6.RIZ ' 02346 16 02412 04988 FLAGS 

16 TFM R.X3+6.6199.9Z 02*6 16 02532 06199 FLAGS 

»7. _ _JZ M _ R 2^2 + Ji ,R l +8Z _ I 023fc0 16 0250 1 04996 FLAGS 

IB TFM RX2+23.R1+88Z 02^82 16 02513 05076 FLAGS 

19 AM *+18.2Z 02394 1.1 02412 00002 FLAGS 

20 RX1 CM R1+2.0.10Z 02406 14 04990 00000 FLAGS 

21 BI CALC.1200Z 02418 46 02586 01200 

22 TF RX3+21 1RXI+6. 1 1 Z 02430 26 02547 02412 FLAGS I ND ADO 
RI >L R2-79.5C0Z _ 02442 36 0610 1 fl.05Q0_. 

2* AM RX3+6.7.9Z 02454 11 02532 00007 FLAGS 

25 AM RX2+U.2Z 02466 1 1 02501 00002 FLAGS 

26 AM RX2+23.2Z 02478 11 02513 00002 FLAGS 

27 RX? TF *+30.0Z 02490 26 02520 00000 

28 TF *+35.0Z 02502 26 02537 00000 

29 SF R2Z _ _ 0251 4 32 06180 00000 

30 RX3 TF 6206. R2Z 02526 26 06206 06180 

31 SM *+9,l«10Z 02538 12 02547 00001 FLAGS 

32 BNI RX2-36.1200Z 02550 47 02454 01200 

33 CM RX1+6,R1+8Z 02562 14 02412 04996 FLAGS 

34 BNI RX1-12.1200Z 02574 47 02394 01200 

35 CALC _ TFM CX 1 +6 . 6479 > 9Z _ Q2586 1_6 0270J3 06479 FLAGS 

36 TFM CX1 + 1 1 .6199.9Z " -■ ■ - ■ — - 02598 1 6 02705 06199 FLAGS 

37 TFM CX1+30.6879,8Z 02610 16 02724 06879 FLAGS 

38 TFM CX1+47.0.10Z 02622 16 02741 00000 FLAGS 

39 TFM CX2+18.7519Z 02634 16 02832 07519 FLAGS 

40 AM CX1+47.1.10Z 02646 11 02741 00001 FLAGS 

41 AM CX1+11.7.9Z 02658 11 02705 00007 FLAGS 

42 AM CX1+30. 16. BZ 02670 11 02724 00016 FLAGS 

43 AM CX1+6.10.9Z 02682 11 02700 00010 FLAGS 

44 CXI A 6489.6206Z 02694 21 06489 06206 

45 M *-l,*-1.611Z 02706 23 0270b 02705 FLAGS 1 ND ADD 

46 A 6B95i99Z 0271S 21 06895 00099 

47 CM R1.1.10Z 02730 14 04988 00001 FLAGS 
4B Bl CX3.1200Z 2742 "46 02874 01 200 " 

49 TF CX2+1 1.CX1+1 1Z 02754 26 02825 02705 

50 TF C.X2+35.CX1+47Z 02766 26 02849 0274 1 

51 AM CX2+11.7.9Z 02778 11 02825 00007 FLAGS 

52 AM CX2+18.16Z . 02790 11 02832 00016 FLAGS 

53 AM CX2 + 35.1.10Z 02802 11 02849 0000 1 FLAGS 

54 CX2' M CX 1 + 1 1 , 62 1 3. 6Z '" 02814 23 02705 06213 FLAGS I ND ADD 

55 A 7519. 99Z 02826 21 07519 00099 



Rl .2. 10Z 
CX2-36. 1 200Z 
CALC+60Z 
NtlZ 

READ 1 900Z 
N'M'i . NZ 
NM1 , 1 Z 
FLOAT .NM 1 Z 
PUNCH 1 . 1 OOZ 
. 1 02Z 
T I TLE1 . 1 007 
. 1 02Z 

MSD1 + 1 1 . 64 79 1 
MSD2+6.6879Z 
PMD1 < • 1 OZ 
MSD1+1 1 • 10. 9Z 
MSD2+6. 1 6Z 
PMDl'. 1 . 1 OZ 
9 7 .6489Z 
9 1 .NZ 
* .94Z 
94.5. 10Z 
93 . *-24Z 
PMD2 . 93Z 
MSD1 +11. 
79, OZ 
84, NZ 
6895, 94Z 
FLOAT , *- 
99, NM 1 Z 



FSQR PMD3.99Z 



BNF 

TDM 

CF 

TNF 

TNF 

TNF 



SOI +1 1 .61 1Z 



102 
1 03 

104 MSD3 



*+23.PM02-7Z 

*+60 ,Z 

*- 1 , PMD2- 1 Z 

* + 36 . 1 200Z 

*-25. 1Z 

*-48Z 

*-49.0.6Z 

PMDA 1 +48 . Bl GB+4 8Z 

*+36 . PMD2Z 

PMDA 1 +23, 2Z 

PMD2Z 

PMDA 1 +2 , PMDl Z 
~P'MDAV+22,PMD2Z 
PMDA 1 +48,PMD3Z 
PUNCH2.100Z 
. 102Z 
PMDA 1 . 1 OOZ 
PMDl iRIZ 
M3DT-36.T2T3TT2 
PUNCH3. 100Z 



02838 
02850 
02862 
02874 
02886 
02898 
02910 
02922 
02934 
02946 
02958 
02970 
02982 
02994 
3006 
030 1 8 
3030 
3042 
' 03054 
3066 
03078 
03090 
03102 
3 114 
03126 
03138 
03150 
03 162 
03 174 
03186 
03198 
0321 
03222 
03234 
03246 
03258 
03270 
3282 
03294 
03306 
0331 8 
03330 
03342 
03354 
3366 
03378 
03390 
" 6 340 2" 
341 4 
03426 
03438 
03450 
03462 



14 04988 
47 02778 
49 02646 

1 1 04976 
47 02346 
26 04986 

12 04986 
1 7 04062 
46 04446 
34 00000 
39 04887 
34 00000 
16 0306b 
16 03 168 
16 Obi -38 

I 1 C3065 

II 03 1 68 
11 05158 
2b 00097 
2y 00091 
7 1 030 78 
1 1 00094 
71 00093 

26 05166 
23 03065 
16 00079 
29 00084 
22 06895 

27 04062 
16 00469 
26 01260 
49 01422 
16 05704 
16 05692 
49 05178 
\b 0328 1 

43 03330 

14 03281 
46 03330 
1 1 0328 1 
49 03270 

32 03281 
26 04723 

44 03390 

15 04698 

33 05166 
73 04677 
73 04697 
73 04723 
46 04506 

34 00000 
39 04675 
24 05158 



03474 
03486 



47 030 18 
46 04554 



000 2 
1 200 
00000 
OO'OO 1 
00900 
04976' 

0000 1 
C4986 

001 00 
00102 
00100 

102 
064 79 
06879 
00000 

0001 
0001 6 
0000 1 
06489 
049 76 
00094 
0000b 
03078 
0093 
03065 
OOOOC 
0497o 
00094 

03 1 68 
0322 1 
0099 
04986 
03257 
O b 1 7 a 
00099 
051 59 
00000 
051 65 
1 200 

0000 1 
00000 
00000 

04 775 
051 66 
00002 
00000 
05158 

"05166 
051 76 

001 00 
00102 
001 00 
04988 

"6T200" 
00100 



FLAGS 
FLAGS 



FLAGS 
FLAGS 
FLAGS 
FLAGS 
FLAGS 
FLAGS 



F L AGS 
FLAGS 



F LAGS 

FLA.GS 
FLAGS 
FLAGS 
FLAGS 
FLAGS 

FLAGS 

FLAGS 

FLAGS IND ADD 



J_07_ 
108 ' 

I 09 

II 

1 1 1 

112 
1 1 3 
1 1 4 
115 
116 

1 1 7 

118 CR1 
1 1 9 
120 
.12 1 
1 22 

123 

1 24 
125 
1 26 

127 CR2 
1 27 

127 

128 
1 28 
128 

129 CR3 
1 30 

131 

1 32" CR4 
1 33 
134 
1 34 
1 34 

1 35 

__ 

1 37 
1 38 
1 39 
1 40 
141 

T42 

1 43 

1 44 

1 45 
.146 

1 47 
S'4fi " r L A ' 

1 49 

1 50 

1 5 1 

1 52 

1 53 FA 

1=4 

1 55 



TFM 
TFM 
TFM 



0. 102Z 

C . 102Z " ' 

TITLE2 .J00Z 
. 1 C2Z 

PMD1 .0.10Z 

CR2+23.6879Z 
CR4+6.7519Z 
CR3+6.64 79.9Z 
PMD1 » 1 « 1 OZ 
CR2 + 23. 16Z 

CR3+ 6 . 10. 9Z 

PMDli'RlZ 
INIT. 1200Z 
CR2+35.CR2+23Z 
CR3+1 1 .CR3+6Z 
PMD2.PMD1 z 
CR2+35.1 6Z 
CR3+1 1 « 1 0.9Z~ 
CR4+6« 16Z 
PMD2. 1 , 1 OZ 
6895.691 1 Z 



FSOR PMD3.99Z 



6489.6499Z 
79, OZ 
84.NZ 
7535 t94Z 
FLOAT, *-6Z 
99.PMD3Z 



8NF 

TDM 

CF 

TNF 

TNF 

TNF 



TFM 
TFM 



PMDA 1 +48 . B 1 GB+4BZ 
- * + 36*.99Z 
PMDA1+35.2Z 
99Z 

PMDA1+34.99Z 
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